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(57) Abstract 

The invention relates to a method for designing a technical system. According to said method, a model of said system is produced 
in a descriptive format An interference is modelled using a random variable. The system, including the interference, is modelled in the 
form of a stochastic differential equation. This stochastic differential equation is solved numerically and the solution obtained is used to 
design the technical system. 



(57) Zusammenfassung 

Es wird ein Verfahren zum Entwurf eines technischen Systems angegeben, bei dem das System in Form einer Beschreibungsform 
modelliert ist. Eine Stoning wird anhand einer ZufallsgroBe modelliert. Das System einschlieBlich der Stoning wird als stochastische 
Differential gleichung modelliert. Die stochastische Differential gleichung wird numerisch gelost und die ermittelte Losung zum Entwurf des 
technischen Systems eingesetzt. 
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Beschreibung 

Verfahren und Anordnung zum Entwurf eines technischen Systems 

5 Die Erfindung betrifft ein Verfahren und eine Anordnung zum 
Entwurf eines technischen Systems, 

Fur den Entwurf eines komplexen technischen Systems ist es 
notwendig, unter einer Menge zulassiger Arbeit spunk te (Be- 

10 triebspunkte, Auslegungsparameter) mindestens einen Ar- 

beitspunkt zu bestimmen, der einen effektiven Ablauf des Sy- 
stems ermoglicht. Dabei sind oftmals Nebenbedingungen zu be- 
rticksichtigen, die die Menge der Arbeitspunkte einschranken 
und dementsprechend bei der Suche nach dem Arbeitspunkt mit- 

15 zuberucksichtigen sind. 



In der Praxis wird das technische System haufig anhand der 
Erfahrung und des Wissens eines (oder mehrerer) Experten ent- 
worfen. Dabei ist es von Nachteil, daJ5 bei Oberschreitung ei- 
20 ner gewissen Komplexitat fur das technische System die 

Fehlerwahrscheinlichkeit beim manuellen Entwurf ansteigt. 

Aus [1] sind Brownsche Bewegungen auf Mannigf altigkeiten und 
insbesondere deren Darstellung mit Hilfe von Projektionsope- 
25 ratoren und Stratonovich-Integralen bekannt. 

Das Konzept einer Stratonovich-Gleichung (auch: "Fisk- 
Stratonovich-Gleichung") ist aus [3] bekannt. Dabei bezeich- 
net die Stratonovich-Gleichung einen Typ stochastischer Dif- 
30 ferentialgleichungen mit speziellen (siehe [3]) Eigenschaf- 
ten. 

Nun unterliegt gerade besagter Entwurf oftmals einem physika- 
lischen EinfluB, der als Storung identif izierbar ist. Diese 
35 Storung wirkt sich u.U. signifikant auf den Entwurf aus, so 
dafi auf der Suche nach einer akzeptablen Losung fur den Ent- 
wurf die Storung nicht unberucksichtigt bleiben kann. 
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Die Aufgabe der Erfindung besteht darin, einen Entwurf eines 
technischen Systems zu ermoglichen, das einer Storung unter- 
liegt . 

Diese Aufgabe wird gemafl den Merkmalen der unabhangigen Pa- 
tentanspriiche gelost. Weiterbildungen der Erfindung ergeben 
sich auch aus den abhangigen Anspriichen. 

Zur Losung der Aufgabe wird ein Verfahren zum Entwurf eines 
technischen Systems angegeben, bei dem das System in Form ei- 
ner Beschreibungsform modelliert ist. Eine Storung wird an- 
hand einer Zuf allsgroiie modelliert. Das System einschlieBlich 
der Storung wird als stochastische Dif f erentialgleichung mo- 
delliert. Die stochastische Dif f erentialgleichung wird nume- 
risch gelost und die ermittelte Losung zum Entwurf des tech- 
nischen Systems eingesetzt. 

Eine Weiterbildung besteht darin, daiJ die Beschreibungsform 
eine Dif f erentialgleichung ist. 

Eine Umsetzung des Ergebnisses aus Modellierung bzw. Simula- 
tion erfolgt vorzugsweise, indem anhand der Losung der sto- 
chastischen Dif f erentialgleichung mindestens ein akzeptabler, 
wenn moglich sogar ein effizienter Arbeitspunkt des techni- 
schen Systems ermittelt wird. Als effizienter Arbeitspunkt 
wird dabei ein zulassiger Arbeitspunkt des Systems bezeich- 
net, der insoweit optimal ist, als eine Veranderung dieses 
Punktes keine Verbesserung des Systems mehr ermoglicht. Hier- 
bei sei darauf hingewiesen, dafi der Arbeitspunkt mehrerer Pa- 
rameter umfassen kann, deren jeder eine einstellbare Grofie 
sein kann. 

Ein Vorteil besteht darin, daJ5 durch die stochastische Diffe- 
rentialgleichung die Storung modellierbar wird und damit im 
Entwurf des technischen Systems mitberucksichtigt werden 
kann. Gerade bei realen technischen Systemen, bspw. in der 
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Schaltkreissimulation, spielen reale Storungen eine bedeuten- 
de Rolle: So kann sich aufgrund der Storung das Verhalten ei- 
ner Schaltung grundlegend verandern bzw. die Schaltung u.U. 
sogar nicht mehr funktionieren. Durch eine explizite Beruck- 
5 sichtigung der Storung ist es moglich, einen fur die Storung 
angepaBten Entwurf des technischen Systems vorzunehmen. 

Insbesondere ist es eine Ausgestaltung, daB die stochastische 
Dif ferentialgleichung eine Stratonovich-Gleichung ist. Da- 
10 durch wird es moglich, eine ggf. in Form einer (deterministic 
schen) Dif ferentialgleichung vorliegende Modellierung des Sy- 
stems im Hinblick auf vorhandene Storungen derart anzupassen, 
daB ein Entwurf des Systems einschlieBlich Storung moglich 
ist. 

15 

Als eine Weiterbildung ist die Stratonovich-Gleichung derart 
bestimmt, daB deren Losung in einem Def initionsbereich der 
Stratonovich-Gleichung (d.h. in einer gekriimmten Flache, ei- 
ner sogenannten 'Mannigf altigkeit 1 ) verlauft. 

20 

Eine andere Weiterbildung besteht darin, daB die Stratono- 
vich-Gleichung mittels eines Pradiktor-Korrektor-Verf ahrens 
numerisch gelost wird. Ein Pradiktor ist allgemein beschreib- 
bar als eine VorhersagegroBe, die gegeniiber der zu vorhersa- 
25 genden GroBe einen Fehler aufweist, welcher mittels des Kor- 
rektors weitgehend ausgeglichen wird. 

Insbesondere kann der Pradiktor anhand der Trapezregel be- 
stimmt werden. Dabei ist die Trapezregel fur Riemann- 

30 Integrale beispielsweise aus [4] bekannt. Die Trapezregel fur 
Stratonovich-Integrale folgt aus der Definition des Stratono- 
vich-Integrals selbst. So ist das Stratonovich-Integral als 
ein Grenzwert der Trapezsummen definiert, die Trapezregel ist 
damit ein Teil der Definition des Stratonovich-Integrals . Es 

35 wird mittels Trapezsummen eine implizite Gleichung fur den 

Pradiktor aufgestellt, die anhand eines Newton-Verf ahrens ap- 
proximate gelost wird. 
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Eine Ausgestaltung besteht darin, dafi zusatzlich eine 
Schrittweitensteuerung im Verlauf der numerischen Berechnung 
erfolgt. Dies geschieht insbesondere durch eine Schrittwei- 
tenverkleinerung, falls ein zu grower Approximationsf ehler 
vermutet wird (Wert fur den vermuteten Approximationsf ehler 
iiberschreitet eine vorgegebene Schranke) , oder durch eine 
SchrittweitenvergroBerung, falls ein (sehr) kleiner Approxi- 
mationsf ehler vermutet wird (Wert fur den vermuteten Approxi- 
mationsf ehler liegt unterhalb einer vorgegebenen Schranke) . 

Eine andere Ausgestaltung besteht darin, daii die Stratono- 
vich-Gleichung insbesondere die folgende Form aufweist: 

t 

x t = x a - Jp(x s ) ■ F(X S ) . ds + Z (1) 
a 

mit t > a und X a = x Q 
wobei 



a 


einen Startzeitpunkt , 


*0 


einen zulassigen Startparametervektor, 


Xt 


einen Losungsparametervektor zum Zeitpunkt t, 


P(X S ) 


eine Projektionsmatrix im Parametervektor X s , 


F(X S ) 


eine Drif trichtung im Parametervektor X s , 


z 


eine Zuf allsgrbiie 



bezeichnen. 

In einer zusatzliche Ausgestaltung ist die Zuf allsgrofie Z be- 
stimmt durch 

t 

Z= E • J P(x s ) o d B s (3)/ 
a 



wobei 

8 



einen Skalierungsf aktor und 
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{b s : s £ a} eine n-dimensionale Brownsche Bewegung 
bezeichnen. 

Liegt zusatzlich das technische System in Form von Gleichun- 
gen und Ungleichungen vor, die Punkte des technischen Systems 
eingrenzen, so wird insbesondere zunachst das technische Sy- 
stem mittels Einfiihrung von Schlupfvariablen so umf ormuliert , 
daft es eine Mannigf altigkeit beschreibt. In der obigen Diffe- 
rentialgleichung (4) werden nur noch diejenigen Anteile der 
Integranden beriicksichtigt, die tangential zur Mannigf altig- 
keit verlaufen. Dies geschieht durch Multiplikation mit Pro- 
jektionsmatrizen. Es ergibt sich die besondere stochastische 
Dif ferentialgleichung nach Gleichung (1) (insbesondere in Zu- 
sammenschau mit den Gleichungen (2) und (3)), die vom Typ ei- 
ner Stratonovich-Dif ferentialgleichung ist. 

Eine zusatzliche Weiterbildung besteht darin, daft eine das 
technische System beschreibende Ungleichung h-j^x) < 0 mittels 
einer Schlupf variable umgeformt wird zu einer Gleichung. Die- 
se Umformung hat insbesondere die Form 



^i(x, s) = hi(x) + g(s) 



(5), 



wobei 



g(s) 



eine Funktion mit einem Wertebereich aller reel- 
len Zahlen groBer gleich Null ( Rq ) / 
die Schlupfvariable und 
eine Nebenbedingung h^(x) < 0 



s 



hi(x) 



bezeichnen. 



Ferner kann die Funktion g(s) bestimmt sein zu: 



g(s) = s 2 



(6) . 



Die Losung der Gleichung (1) stellt einen stochastischen Pro- 
zeJi 
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{x t :t > a} (7) 

dar. Es wird eine Folge zulassiger Parametervektoren 

x ty x t 2 ' x t3/ •• • ( a < ti < t 2 < t 3 . . .) (8) 

berechnet, die in Naherung Punkte eines Pfades des stochasti- 

schen Prozesses gemaB Gleichung (7) zu den Zeitpunkten 

ti, t2* t3, • . . darstellen. Ein Iterationsschritt kommt dabei 

auf die folgende Art zustande: Die Parametervektoren 

x ty x t 2 ' x t 3 ' • • • f x ti seien bereits berechnet, x ti dient dann 

als Startpunkt fur die Stratonovich-Gleichung (1), d.h. 
a:= ti . Ein Parametervektor x ti + 1 wird als numerische Ap- 
proximation eines Punktes des Pfades der Losung gemali Glei- 
chung (7) , der bereits die Punkte x tl , x t2 , x t3 , . . . , x ti ent- 

halt, berechnet. Die Berechnung erfolgt zum Zeitpunkt ti +1 , 
d.h. t = ti+i, die Schrittweite des Iterationsschritts er- 
gibt sich zu t± + i - t^ . 

Die Approximation der Stratonovich-Gleichung (1) erfolgt fiir 
den (i+l)-ten Iterationsschritt nach der folgenden Vorschrift 
(mit Hilfe der Trapezregel) : 

X ti+l - x ti 

-H±^^ t J. F (x ti ) + p(x ti+1 ). F (x ti + 1 )) 

+ f ^ti)+^t i+1 )-tt i + 1 -B t J 

(9) 

Gleichung (9) ist eine implizite Gleichung fur die Zufallsva- 
riable X ti + 1 . Mit Hilfe des Newton-Verf ahrens kann daraus 
eine Approximation yt i + 1 fiir eine Realisierung der Zufalls- 
variablen X^. . bestimmt werden. 
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Der Punkt Yti +1 ist u - u - nicht zulassig* Anhand des Korrek- 

tors wird nun mit Hilfe eines Newtonverf ahrens ein Punkt 
Xt i+1 in der N * he des p unktes yt^ +1 ermittelt, welcher zu- 
lassig ist und als Startwert fur einen nachfolgenden Iterati- 
5 onsschritt dient. 

Eine andere Ausgestaltung besteht darin, dafi der Entwurf ein 
Neuentwurf, eine Anpassung oder eine Steuerung des techni- 
schen Systems darstellt. Alternativ kann im Rahmen des Ent- 
10 wurfs eine Modeliierung und/oder Simulation des technisches 
Systems erf oigen. 

Ein technisches System kann eine Anlage der Verf ahrenstechnik 
oder ein sonstiges System sein, das im Hinblick auf unter- 
schiedliche Parameter auszulegen oder einzustellen ist. Ins- 
besondere konnen die Parameter des Parametervektors Ausle- 
gungsparameter oder Betriebsparameter des technischen Systems 
sein. Betriebsparameter kennzeichnen mogliche einstellbare 
Groflen, wohingegen Auslegungsparameter insbesondere physika- 
lische Abmessungen des technischen Systems beschreiben und 
wahrend des Betriebs zumeist nur mit hohem Aufwand angepaflt 
oder verandert werden konnen. 

Wahlweise kann mit dem beschriebenen Verfahren ein Neuentwurf 
25 des technischen Systems oder eine Anpassung eines bereits 

existierenden technischen Systems erfolgen. In beiden Fallen 
handelt es sich urn einen Entwurf (einmal als Neugenerierung 
und einmal als Anpassung) im Sinne der vorliegenden Ausftih- 
rungen . 

30 

Im Rahmen einer weiteren Ausgestaltung wird das technische 
System anhand der ermittelten Parameter realisiert bzw. ein- 
gestellt. Hierbei ist es von Vorteil, daii die Parameter in 
einen Parametervektor , der mittels der Erfindung bestimmt 
35 wurde, einen stabilen Betriebspunkt kennzeichnen und die Ein- 
stellung des Systems auf diesen Betriebspunkt einen dauerhaft 
sicheren Betrieb des Systems/ der Anlage gewahrleistet . 



15 



20 
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Auch wird die Losung anhand einer Anordnung zum Entwurf eines 
technischen Systems gelost, die eine Prozessoreinheit auf- 
weist, die derart eingerichtet ist, dafl 

a) das System in Form einer Beschreibungsf orm modelliert 
ist; 

b) eine Storung des Systems anhand einer Zuf allsgroBe mo- 
dellierbar ist; 

c) das System einschlieBlich der Storung als eine stocha- 
stische Dif f erentialgleichung modellierbar ist; 

d) die stochastische Dif f erentialgleichung numerisch 16s- 
bar ist; 

e) die numerische Losung zum Entwurf des technischen Sy- 
stems erfolgt. 

Diese Anordnung ist insbesondere geeignet zur Durchfiihrung 
des erf indungsgemafien Verfahrens oder einer seiner vorstehend 
erlauterten Weiterbildungen. 

Ausfuhrungsbeispiele der Erfindung werden nachfolgend anhand 
der Zeichnung dargestellt und erlautert. 

Es zeigen 

Fig.l einen Ablauf eines Verfahrens zum Entwurf eines tech- 
nischen Systems; 

Fig. 2 eine Skizze, die eine numerische Losung einer Strato- 
novich-Diff erentialgleichung veranschauiicht; 

Fig. 3 eine Skizze, die eine Projektionsmethode veranschau- 
iicht; 



Fig. 4 eine Prozessoreinheit. 
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Fig.l zeigt ein Blockdiagramm, das den Ablauf eines Verfah- 
rens zum Entwurf eines technischen Systems veranschaulicht. 
In einem Schritt 101 ist eine Beschreibungsf orm, bevorzugt 
als deterministische Dif f erentialgleichung bzw. als determi- 
nistisches Dif f erentialgleichungssystem, eines technischen 
Systems vorgegeben. In einem Schritt 102 wird eine Storung, 
die auf das technische System einwirkt, berucksichtigt, indem 
eine Stratonovich-Gleichung formuliert und in Schritt 103 ge- 
lost wird. Die Losung der Stratonovich-Gleichung fiihrt gemaJi 
Schritt 104 zu einem Simulationspunkt fur eine dem techni- 
schen System zuzuordnenden Grofle unter Berlicksichtigung der 
Storung. Eine sich ggf. mit dem Simulationspunkt anschliefien- 
de Optimierung fUhrt zu einem Arbeitspunkt (insbesondere ei- 
nem effizienten Arbeitspunkt), der z.B. fur die Umsetzung des 
technischen Systems eingesetzt wird. 

Fig. 2 veranschaulicht das Pradiktor-Korrektor-Verf ahren. Aus- 
gehend von einem Punkt / der sich im zulassigen Wertebe- 

reich L befindet, wird ein Punkt yt i + 1 ermittelt (Anwendung 
des Newton-Verfahrens auf Gleichung (9)), der nicht mehr im 
zulassigen Bereich L liegt (Pradiktorschritt 201). Urn den 
Punkt Yti +1 wieder in den zulassigen Bereich zu uberfuhren 

wird ein Korrektorschritt 202 bestimmt, der zu dem Punkt 
x ti+l fuhrt - 

Das Zusammenspiel von Pradiktor und Korrektor wird iterativ 
mit gleichbleibender, abnehmender oder zunehmender Schritt- 
weite wiederholt. Hierbei sei bemerkt, daft die Linie 201 
nicht notwendigerweise tangential verlauft. 

Fig. 3 veranschaulicht die Projektionsmethode. Wiederum ist 
der zulassiger Bereich L angegeben, auf dem sich der Punkt 
x ti befindet. Die Linie 203 kennzeichnet eine Drif trichtung 

F(X t ), eine Linie 204 deutet eine projizierte Drif trichtung 
P(x t ) • F(x t ) an. 

In Fig. 4 ist eine Prozessoreinheit PRZE dargestellt. Die Pro- 
zessoreinheit PRZE umfaflt einen Prozessor CPU, einen Speicher 
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SPE und eine Input/Output-Schnittstelle IOS, die iiber ein In- 
terface IFC auf unterschiedliche Art und Weise genutzt wird: 
t)ber eine Graf ikschnittstelle wird eine Ausgabe auf einem Mo- 
nitor MON sichtbar und/oder auf einem Drucker PRT ausgegeben. 
5 Eine Eingabe erfolgt uber eine Maus MAS oder eine Tastatur 
TAST. Auch verfugt die Prozessoreinheit PRZE tiber einen Da- 
tenbus BUS, der die Verbindung von einem Speicher MEM, dem 
Prozessor CPU und der Input/Output-Schnittstelle IOS gewahr- 
leistet. Weiterhin sind an den Datenbus BUS zusatzliche Kom- 
10 ponenten anschlieflbar, z*B. zusatzlicher Speicher, Datenspei- 
cher (Festplatte) oder Scanner. 
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Patent anspriiche 

1. Verfahren zuiti Entwurf eines technischen Systems, 

a) bei dem das System in Form einer Beschreibungsf orm mo- 
delliert ist; 

b) bei dem eine Storung des Systems anhand einer Zufalls- 
grofte modelliert wird; 

c) bei dem das System einschliefilich der Storung als eine 
stochastische Dif f erentialgleichung modelliert wird; 

d) bei dem die stochastische Dif f erentialgleichung nume- 
risch gelost wird; 

e) bei dem die numerische Losung zum Entwurf des techni- 
schen Systems eingesetzt wird. 

2. Verfahren nach Anspruch 1, 

bei dem die Beschreibungsf orm eine Dif f erentialgleichung 
ist . 

3. Verfahren nach einem der vorhergehenden Anspriiche, 

bei dem anhand der Losung der stochastischen Differenti- 
algleichung mindestens ein effizienter Arbeitspunkt des 
technischen Systems ermittelt wird. 

4. Verfahren nach einem der vorhergehenden Anspriiche, 
bei dem die stochastische Dif f erentialgleichung eine 
Stratonovich-Gleichung ist. 

5. Verfahren nach Anspruch 4, 

bei dem die Stratonovich-Gleichung derart bestimmt wird, 
daft deren Losung in einem Def initionsbereich der Strato- 
novich-Gleichung verlauft. 

6. Verfahren nach Anspruch 4 oder 5, 

bei dem die Stratonovich-Gleichung mittels eines Pradik- 
tor-Korrektor-Verf ahrens numerisch gelost wird. 
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7. Verfahren nach Anspruch 6, 

bei dem ein Pradiktor des Pradiktor-Korrektor-Verf ahrens 
anhand der Trapezregel bestimmt wird. 

5 8. Verfahren nach Anspruch 6 oder 7, 

bei dem ein Korrektor des Pradiktor-Korrektor-Verf ahrens 
mittels Newton-Verf ahren bestimmt wird. 

9. Verfahren nach einem der vorhergehenden Anspriiche, 

10 bei dem zusatzlich eine Schrittweitensteuerung erfolgt. 

10. Verfahren nach Anspruch 9, 

bei dem die Schrittweitensteuerung derart erfolgt, 

a) daJi eine Schrittweitenverkleinerung durchgefuhrt wird, 
15 falls ein hinreichend groBer Approximationsf ehler er- 

mittelt wird oder 

b) dafl eine SchrittweitenvergroBerung durchgefuhrt wird, 
falls ein hinreichend kleiner Approximationsf ehler er- 
mittelt wird. 

20 

11. Verfahren nach einem der Anspriiche 4 bis 9, 

bei dem die Stratonovich-Gleichung die folgende Form auf- 
weist : 

t 

25 X t = X a - J P(X S ) • F(x s ) ■ ds + Z mit t > a und X a = x 0 

a 

wobei 

a einen Startzeitpunkt (a > 0), 

xo einen zulassigen Startparametervektor, 

30 Xt einen Losungsparametervektor zum Zeitpunkt t, 

P(X S ) eine Projektionsmatrix im Parametervektor X s , 

F(X S ) eine Drif trichtung, 

Z eine Zuf allsgroiie 

bezeichnen. 



35 
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12. Verfahren nach Anspruch 11, 

bei dem die Zufallsgrofie bestimmt ist durch 

t 

Z = 8 • J P(X S ) o dB s 
a 

wobei 

e einen Skalierungsf aktor und 

{b s : s > a} eine n-dimensionale Brownsche Bewegung 

bezeichnen. 

13. Verfahren nach einem der vorhergehenden Anspruche, 

bei dem eine Nebenbedingung, die in Form einer Unglei- 
chung vorliegt, mittels einer Schlupfvariable s in eine 
Gleichung umgef ormt wird. 

14. Verfahren nach Anspruch 13, 

bei dem die Nebenbedingung hi (x) , die in Form der Unglei- 
chung hi(x)<0 vorliegt, formuliert wird zu 

hi(x, s) = hi(x) + g(s), 

wobei 

g(s) eine Funktion mit einem Wertebereich aller reel- 

len Zahlen grofier gleich Null, 
s die Schlupfvariable und 

hi(x) eine Nebenbedingung h^(x) < 0 
bezeichnen. 

15. Verfahren nach Anspruch 14, 

bei dem die Funktion g(s) bestimmt ist zu 



gfe) = s 2 . 
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16. Verfahren nach einem der vorhergehenden Anspruche, 

bei dem der Entwurf ein Neuentwurf , eine Anpassung oder 
eine Steuerung des technischen Systems darstellt. 

17. Anordnung zum Entwurf eines technischen Systems 

mit einer Prozessoreinheit, die derart eingerichtet ist, 
daii 

f ) das System in Form einer Beschreibungsf orm modelliert 
ist; 

g) eine Storung des Systems anhand einer Zuf allsgrofte mo- 
dellierbar ist; 

h) das System einschliefllich der Storung als eine stocha- 
stische Dif f erentialgleichung modellierbar ist; 

i) die stochastische Dif f erentialgleichung numerisch 16s- 
bar ist; 

j) die numerische Losung zum Entwurf des technischen Sy- 
stems erfolgt. 
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